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Abstract 

MicroRNAs are endogenous, small non-coding RNAs approximately 18-26 nucleotides in length that regulate target gene 
expression at the post-transcription level. Interferon-y (IFN-y) is a Thl cytokine that is involved in both the innate and 
adaptive immune responses. We previously identified two IFN-y genes in green-spotted puffer fish {Tetraodon nigroviridis). 
To determine whether miRNAs participate in IFN-y-related immune responses, 7. nigroviridis spleen cells were treated with 
recombinant IFN-y isoforms, and a Solexa high-throughput sequencing method was used to identify miRNAs. In total, 1,556, 
1,538 and 1,573 miRNAs were found in the three samples, and differentially expressed miRNAs were determined. In total, 
398 miRNAs were differentially expressed after rIFN-yl treatment, and 438 miRNAs were differentially expressed after rlFN- 
y2 treatment; additionally, 403 miRNAs were differentially expressed between the treatment groups. Ten differentially 
expressed miRNAs were chosen for validation using qRT-PCR. Target genes for the differentially expressed miRNAs were 
predicted, and GO and KEGG analyses were performed. This study provides basic knowledge regarding fish IFN-y-induced 
miRNAs and offers clues for further studies into the mechanisms underlying fish IFN-y-mediated immune responses. 
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introduction 

MicroRNAs (miRNAs) are small non-coding RNAs that 
regulate target messenger RNA (mRNA) expression by binding 
to mRNA 3' untranslated regions (3' UTR), resulting in mRNA 
cleavage or translational repression through the RNA-induced 
silencing complex (RISC) [1]. MiRNAs play an important role in 
many biological processes, such as development [2,3], cell 
proliferation, differentiation and death [4,5,6,7], tumorigenesis 
[8,9], and inflammation [10,11], in addition to several diseases 
[12,13]. In recent years, with the development of high-throughput 
sequencing technologies, miRNAs have been increasingly identi- 
fied in teleosts such as Japanese flounder (Paralichthys olivaceus) [14], 
Atlantic salmon [Salmo salar) [15], Atlantic halibut [Hippoglossus 
hippoglossus) [16,17], and channel catfish {Ictalums punctatus) [18]. 
Several miRNA expression studies have been performed in 
zebrafish, including miRNA profiling during development [19], 
miRNA profiling in response to LPS and bacterial infection [20], 
and miRNA profiling induced by PFOS [21]. 

Interferon (IFN)-y, the only type II IFN [22], regulates immune 
processes, including antigen presentation, antiviral response, 
macrophage activation, and apoptosis [23], by binding to the 
IFN-y receptor (IFNGR) [24]. IFN-y ligand-receptor binding 
results in activation of Janus tyrosine kinase Jak)l [25], followed 



by phosphorylation of the signal transducer and activator of 
transcription (Stat)l. The Stat 1 -Stat 1 homodimer binds to GAS 
(interferon-gamma activated sequence) elements in the promoters 
of IFN-y-regulated genes and leads to gene regulation [26]. Stat3 
and Stat5 can also be activated by IFN-y in certain cell types 
[27,28]. In teleosts, IFN-y was first identified in Fugu [Takifugu 
rubripes) in 2004 [29]. Fish IFN-y induces typical IFN-y-induced 
genes, including the major histocompatibility complex (MHC) 
class II and Statl, and also primes macrophages and enhances 
respiratory bursting, suggesting that IFN-y has similar functions in 
fish and mammals [30,31,32]. In contrast to the single copy of the 
IFN-y gene that exists in mammalian genomes, two copies of the 
IFN-y gene were found in the zebrafish and fugu genomes in 2006; 
these genes were named IFN-y 1 and IFN-y 2 [33]. The IFN-y 2 
gene has proven to be an orthologxie of mammalian IFN-y, 
whereas IFN-y 1 shares low identity with other species and appears 
to be fish-specific [34]. Two IFN-y genes were also found in the 
channel catfish (/. punctatus) [35], common carp {Cyprinus carpio L.) 
[36], goldfish {Carassius aurutus L.) [37], and ginbuna crucian carp 
[Carassius auratus langsdorfii) [38]. Both isoforms are members of the 
type II interferon family and contain IFN-y signature motifs 
[30,33], but the fish-specific IFN-yl does not contain the C- 
terminal nuclear translocation signal (NLS) motif required for 
IFN-y activity [35,39]. Previous studies showed that fish IFN-yl 



PLOS ONE I www.plosone.org 



1 



May 2014 | Volume 9 | Issue 5 | e96336 



Tetraodon Mirome after Two IFN-g Stimulations 



(A) 



Length Distribution 



0.01 0.02 0.03 0.04 0.06 0.10 




102 osa 0.44 0.30 0.08 0.03 Q.Qi 



10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 

Length(nt) 



(B) 



Length Distribution 



50.13 












22. 

1 


9 


11.23 






4.27 

■ Hi ^ 0.6S 0.36 0.20 0.05 0.02 



(0 



10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 

Length(nt) 

Length Distribution 



10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 

Length{nt) 



Figure 1. Read length distribution of tKie clean reads from the (A) Sp-con, (B) Sp-yl, and (C) SP-y2 groups. 

doi:10.1371/joumal.pone.0096336.g001 



and IFN-Y2 exhibited different expression patterns [33,39] and 
distinct preferences for different receptors [38,40], but both 
induced typical IFN-y target genes [41], indicating that the two 
type II IFNs may have separate functions [40,41]. However, the 
functional differentiation of the two isoforms and how they 
differentially regulate downstream genes are still not clear. Because 
the involvement of miRNAs in the regulation of immune response 
has been extensively demonstrated [42,43], it is possible that 



miRNAs might be involved in fme-tuning the functions of the two 
type II IFN isoforms. 

MiRNAs have been well studied in mammalian inflammation 
and tumors. Because of the importance of IFN-y in the immune 
system, increasing attention has been paid to the mlRNA network 
after IFN-y stimulation, and several miRNAs have already been 
reported to act on mammalian IFN-y. The miR-29 family directly 
targets IFN-y [44], and IFN-y-induced Stall can up-regulate miR- 
29 expression [45]; in cancer cells, miR-145 directly targets Statl 
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Figure 2. Small RNA tags matched to the genome by SOAP. The figure shows the expression and genomic distribution of miRNAs identified in 
the (A) Sp-con, (B) Sp-yl, (C) SP-y2 groups. 
doi:1 0.1 371 /journal.pone.0096336.g002 



and c-Myc [46,47], and miR-378 targets the IFN-y receptor 1, 
IFNGRl, to suppress bovine luteal cell apoptosis [48]. However, 
no study has focused on the miRNA expression profile in response 
to teleostean IFN-y stimulation, and comparing the differential 
expression patterns would provide valuable insight into the 
separate functions of the two type II interferon isoforms in fish. 

The green-spotted puffer fish {Tetraodon nigroviridis) has the 
smallest known vertebrate genome. Litde information regarding T. 
nigroviridis miRNAs can be found in the miRBase, and until now, 
no studies have evaluated T. nigroviridis miRNAs and T. nigroviridis 
IFN-y. We previously identified two IFN-y genes in the T. 
nigroviridis genome, IFN-yl (GenBank accession NO. KJ524454) 
and IFN-y2 (KJ524455). Recombinant IFN-yl (rIFN-yl) and 
IFN-Y2 (rIFN-y2) from T. nigroviridis were produced in Escherichia 
colt using a pET expression vector, as previously described [49], 
and binding of the recombinant proteins to their receptors was 
confirmed using ligand binding analysis [49]. In this study, 
primary' cells from the T. nigroviridis spleen were treated with rlFN- 
yl and rIFN-y2. An lUumina deep-sequencing method was 
performed, and the differential expression of miRNAs was 
analyzed. These findings will help to determine the roles that 
miRNA play in regulating IFN-y-mediated immune responses. 

Materials and Methods 

Ethics statement 

AU animal experiments were conducted in accordance with the 
guidelines and approval of the Animal Research and Ethics 
Committees of Sun Yat-Sen University. AU efforts were made to 
minimize suffering. 

Experiment animals 

Green-spotted puffer fish were purchased from the local market 
and raised for at least one week prior to the experiment. The fish 
were kept in 28°C circulating water and fed daily. Upon initiation 

of the experiment, the fish were anesthetized with 1 % tricaine 
methanesulfonate (MS-222), and their spleens were removed and 
immediately placed in tissue culture medium (RPMI- 1 640 culture 
medium supplemented with 10% fetal bovine serum, 2 mM L- 
glutamine, and 1% penicillin/streptomycin). 

Sample collection and RNA isolation 

Primary spleen cells from T. nigroviridis were obtained as 
previously described, with minor modifications [50]. Briefly, tissue 
culture medium (TCM) was prepared from RPMI- 1640 medium 
(Life Technologies, Gaithersburg, MD, USA) by adding 10% fetal 
bovine serum (Life Technologies), 2 mM L-glutamine (Sigma, St. 
Louis, MO, USA), and l'!'o penicillin/streptomycin (Sigma). Six 
fish were used in this study; their spleens were removed and 
immediately kept in TCM. The spleens were then ground into 
pieces and filtered through a 70- ^im cell strainer (BD Falcon, 
Bedford, NIA, USA). The obtained cells were washed and 
centrifuged with TCM twice, then dividc-d into thrc-e groups (Sp- 
con, Sp-yl, and Sp-y2). The cells were then plated in a 6-well plate 
and incubated at 27°C under 5% CO2 for 6 hours before rIFN-y 
stimulation. Two IFN-y genes were previously identified in the T. 
nigroviridis genome and were named IFN-yl (GenBank accession 
NO. KJ524454) and IFN-y2 (KJ5244,5,5). Recombinant IFN-yl 
(rIFN-yl) and IFN-y 2 (rIFN-y2) from T. nigroviridis was produced 
as previously described [49]; rlFN-yl and rIFN-y2 proteins at a 



final concentration of 10 ng/mL were added to the culture 
medium, and the control group (Sp-con) was treated with TCM. 
The cells were collected 12 hours after rIFN-y stimulation, and 
Trizol (Life Technologies) reagent was immediately added to 
prevent RNA degradation. Total RNA was extracted using Trizol 
reagent according to the manufacturer's instructions. RNA quality 
was examined using an Agilent 2100 Bioanalyzer. Total RNA 
from the same experimental group was pooled prior to the 
construction of the small RNA library. 

Small RNA library construction and sequencing 

Small RNA fragments (18-30 nt) were isolated from total RNA 
using a Novex 15% TBE-Urea gel (Life Technologies). Then, a 5'- 
adaptor (lUumina, San Diego, CA, USA) was ligated to the small 
RNAs, and the products were purified on a Novex 15% TBE- 
Urea gel to remove any non-ligated adaptors. The 5' ligation 
products (36-50 nt) were then ligated to a 3'-adaptor (Illumina) 
and purified on a Novex 10% TBE-Urea gel (Life Technologies). 
Ligation products with adaptors at both ends (62-75 nt) were then 
reverse-transcribed and PCR-amplified using adaptor primers. 
The amplification products were excised from a 6% TBE-Urea gel 
(Life Technologies). The purified cDNA library was used for 
clustering and sequencing analysis using a HiSeq 2000 lUumina 
sequencer at the Beijing Genomics Institute (BGI), Shenzhen, 
China. 

Data analysis 

Sequencing data were uploaded to the NIH Short Read Archive 
(project accession NO. SRP039387). Raw data from the HiSeq 
sequencing were first filtered and cleaned. The low-quaUty reads, 
including reads with no 3' adapter, reads with no insertion, reads 
with 5' adapter, reads smalk-r than 18 nt, and reads with a polyA 
tail, were rem()\ ed. The length distribution of the clean tags was 
then summarized. The clean reads were then mapped to 
miRBase 19 (http://www.microma.org/) and the T. nigroviridis 
genome to identify the conserved miRNAs. The small RNA reads 
were also aligned to repeat associated RNA, exons and introns of 
mRNA, rRNA, scRNA, snoRNA, snRNA, and tRNA. The small 
RNA tags were then annotated using the software tag2annotation 
developed by BGI. 

Identification of differentially expressed miRNAs and 
target prediction 

The miRNA expression profiles of the control and rlFN-y- 
treated groups were compared to identify differentially expressed 
miRNAs. The expression of the miRNAs in the two samples 

(control and treatment) was normalized to determine the 
expression as transcripts per million (TPM). Then, the fold change 
and P-value were calculated from the normEilized expression, and 
a log2-ratio scatter plot was generated. MiRNAs with similar 
expression patterns were clustered using hierarchical clustering, as 
previously described [51]. 

The targets of the annotated miRNAs were predicted using 
Mireap software developed by BGI. The rules used for target 
prediction were based on previously described suggestions [52,53]. 
KOG analysis was performed to classify the targets by their 
functions. The top three targets with the highest scores among the 
differentially expressed miRNAs were analyzed using GO 
enrichment and KEGG pathway analysis. 
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Figure 3. Expression of the identified miRNAs in tKie three 
different samples using hierarchical clustering, l-iighly expressed 
miRNAs are indicated in red, and miRNAs with low expression are 
indicated in green; grey indicates missing data. The absolute signal 
intensity ranged from -4.0 to +4.0. 
doi:10.1371/journal.pone.0096336.g003 



Reverse transcription and quantitative polymerase chain 
reaction (RT-qPCR) 

T. nigroviridis miRNA expression levels were measured as 
previously described [54]. RNA samples that were used in the 
HiSeq sequencing analysis were reverse-transcribed and amplified 
using the Hairpin-it miRNAs RT-PCR Quantitation kit (Gene- 
Pharma, Shanghai, China). The primers used in the reverse 
transcription and real-time PGR reactions were provided by 
GenePharma (Shanghai, China). Real-time quantitative PGR was 
performed on an ABI 7900HT instrument. The reactions were 
incubated at 95°G for 3 min, followed by 40 cycles of 95°G for 12 
sec and 62°C for 40 sec. All reactions were run in triplicate. T. 
nigrovinAis U6 snRNA was used as a control. Relative miRNA 
expression was calculated using the comparative threshold (2'^^^*^') 
method [48]. The Gt values used in the calculation were the 
means of the triplicates, and the data are presented as means ± 
SEM. Significant differences were analyzed using t-tests, and p< 
0.05 was considered statistically significant. 

Results 

High-throughput sequencing of small RNAs 

To identify differentially expressed miRNAs after rIFN-y 
treatment, spleen cells from T. nigroviridis were divided into three 
groups: the control group (Sp-con), the rlFN-y 1 -treated group (Sp- 
yl), and the rIFN-y2-treated group (Sp-y2). All three groups were 
sequenced separately. In total, 11,899,197 raw reads were 
obtained from the Sp-con group, whereas 12,000,000 raw reads 
were obtained from both the Sp-yl and Sp-y2 groups. After the 
low-quality reads were removed, 11,661,792, 11,821,415, and 
11,705,803 clean reads were obtained from the Sp-con, Sp-yl, 
and Sp-y2 groups, respectively. The length distributions of the 
three groups were analyzed (Figure 1); the reads were mainly 22- 
23 nt in length, similar to previous studies [14,55,56]. The clean 
reads were then mapped to the mature miRNAs of all animals in 
miRBasel9, as only 132 T nigroviridis miRNA records could be 
found in miRBasel9. In total, 1556, 1538, and 1573 mature 
miRNA sequences were identified in the Sp-con, Sp-yl, and Sp-y2 
groups, respectively. All miRNAs were mapped to the T. nigroviridis 
genome (Figure 2). All of the identified miRNA sequences are 
listed in Supplemental Table SI. 

Identification of differentially expressed miRNAs 

The known miRNA reads were compared between the Sp-con/ 
Sp-yl, Sp-con/Sp-y2, and Sp-yl/Sp-y2 groups to identify the 
differentially expressed miRNAs. The log2 ratio and p-value were 
calculated, and miRNAs with a log2 ratio >1 and p<0.05 were 
considered to be up-regulated; miRNAs with a log2 ratio < — 1 
and p<0.05 were considered to be down-regulated. The miRNAs 
dififerentially expressed among different experiment groups were 
summarized using hierarchical clustering (Figure 3) and the log2- 
ratio scatter plot (Figure 4). The above analyses indicate that 191 
miRNAs were significantly up-regulated and that 207 miRNAs 
were significantly down-regialated after rIFN-yl stimulation; after 
rIFN-y2 stimulation, 226 miRNAs were significantly up-regulated 
and 2 1 2 miRNAs were significantly down regxilated. In total, 403 
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Figure 4. Log2-ratio scatter plot of the differentially expressed 
miRNAs between (A) the Sp-con and Sp-yl groups, (B) the Sp- 
con and Sp-Y2 groups, and (C) the Sp-y1 and SP-y2 groups. Each 
point in the figure represents a miRNA. The X and Y axes show the 
expression levels in the two samples. The red points represent miRNAs 
with a ratio >2, the blue points represent miRNAs with a ratio >1/2 and 
£2, and the green points represent miRNAs with a ratio £1/2. 
doi:10.1371/journal.pone.0096336.g004 

miRNAs were differentially expressed between the two groups 
treated with rIFN-yl and rIFN-Y2 (Supplemental Table S2). qRT- 
PCR was used to verify the expression profile of the differentially 
expressed miRNAs. In total, 10 differentially expressed miRNAs 
(tni-miR-106a-3p, tni-miR- 1 24b-3p, tni-miR-132-3p, tni-mlR- 
142-3p, tni-miR-145b, tni-miR-17, tni-miR-223-5p, tni-miR- 
29d-3p, tni-miR-346, and tni-miR-378) were selected for PGR 
validation. The normalized read counts of the 1 0 selected miRNAs 
are shown in Figure 5. Whereas tni-miR- 1 24b-3p and tni-miR-17 
were differentially expressed after both rlFN-yl and rIFN-Y2 
treatment, the other eight miRNAs were differentially expressed 
only after either rIFN-yl or rIFN-Y2 treatment. The real-time 
PGR validation of these 10 miRNAs is shown in Figure 6. AU of 
the 1 0 selected miRNAs showed the same pattern of increase or 
decrease suggested by the sequencing data, although the 
magnitude of the change in tni-miR- 106a-3p, tni-miR-17, and 
tni-miR-346 expression was different between the sequencing 
results and PGR data. 

Target prediction 

The T. nigroviridis genome database was used to predict miRNA 
targets. In total, 165,369 genes were predicted to be targets of 
expressed miRNAs, and these genes were annotated in the KOG 
database and divided into 25 groups (Figure 7). Among these 25 
groups, genes involved in signal transduction mechanisms were the 
most commonly targeted genes (13.83%), followed by genes 
involved in translation, ribosomal structure and biogenesis 
(12.69%) and those involved in the cytoskeleton (11.08%). 

To further investigate the differentially expressed miRNAs, the 
top three targets with the highest scores among all the differentially 
expressed miRNAs were analyzed using GO and KEGG analyses 
(Supplemental Table S3-S6). The results showed that many target 
genes were deeply involved in the IFN-y-mediated immune 
response. For example, SOCS3 (GenBank accession NO. 
GR675748), GALT (CR677385), CRFA22 (AY374494), CISH 
(EF195753), and CRFA28 (CR734177) are involved in the Jak- 
Stat signaling pathway, whereas BPI (CR648786), SYNGR2 
(GR663898), CRFB3 (CR726689, GR733157) are involved in the 
Toll-like receptor signaling pathway. NLRC3 (GR655744), 
HSP90AA (GR646339), and lAP (GR709695) are involved in 
the NOD-like receptor signaling pathway, and MHC I 
(GR635680), IgL-3 (AJ575617), caketiculin (CR639214), and 
GD4-2 (GR733985) are involved in antigen processing and 
presentation. The predicted targets for the differentially expressed 
miRNAs are summarized in Supplemental Table S7. The target 
prediction for the difiTerentially expressed miRNAs and the 
functions of the predicted targets require further investigation. 

Discussion 

MiRNAs have been shown to be involved in the innate and 
acquired immune response [57,58]. Using high-throughput 
sequencing technology, we were able to identify differentially 
expressed miRNAs in T. nigroviridis spleen cells after stimulation 
with recombinant IFN-ys. IFN-y is critically involved in innate 
and adaptive immunity, and is an important activator of 
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Figure 5. Normalized read counts for selected miRNAs. (A): Normalized read counts for miR-124b-3p, miR-142-3p, miR-145b, miR-17, miR-29d- 
3p, and miR-346 following IFNyl treatment, compared to the control group. (B): Normalized read counts of miR-106a-3p, miR-124b-3p, miR-132-3p, 
miR-17, miR-223-5p, and miR-378 following IFNy2 treatment, compared to the control group. 
doi:1 0.1 371/journal.pone.0096336.g005 



macrophages. IFN-y has been studied in teleosts [59], as two 
copies of the IFN-y gene have been found in many fish species 
[35,36,60], thus raising interest in the teleost IFN-y system. T. 
nigroviridis has the smallest known vertebrate genome; however, few 
studies concerning T. nigroviridis miRNAs have been reported [61]. 
Our group previously identified two IFN-y genes in T. nigroviridis, 
and our present study has identified miRNAs involved in the 
response to two IFN-ys. To our knowledge, this is the first study to 
investigate IFN-y-related miRNAs in teleosts. 

Using Solexa high-throughput sequencing, we identified 1556, 
1538, and 1573 mature miRNA sequences in the SP-con, SP-yl 
and SP-y2 samples, respectively. All of the identified miRNAs 
were mapped to miRBasel9 to identify the conserved miRNAs. 
Because only 1 32 miRNA records for T. nigroviridis were found in 
miRBasel9, we mapped the sequencing reads to miRNAs of all 
animals in miRBase. We identified many T. nigroviridis miRNAs 
tiiat were not recorded in miRBase, such as miR-764-5p (5'- 
GGUGCCCGCAUCCUCCUCCA-3'), miR-346 (5'- TGTTG- 
CCCGCATCCTCCAC-3'), and miR-146a-5p (5'-TGAGAACT- 
GAATTCCATCGCTGGTT-3'). AU of the identified T. nigrovir- 
idis miRNAs are listed in Supplemental Table SI. Among these 
miRNAs, 397 miRNAs were differentially expressed between the 
control and IFN-y 1 groups, and 438 miRNAs were differentially 
expressed between the control and IFN-y2 groups; additionally, 
403 miRNAs were differentially expressed between the two IFN-y 
treatment groups (Supplemental Table S2). Different miRNA 
expression patterns were observed after rIFN-yl and rIFN-y2 
treatment, suggesting that both IFN-y 1 and IFN-y 2 may be 
functional and may regulate the T. nigroviridis immune system. 
Next, we verified the differential expression data using qRT-PCR. 
Ten miRNAs were selected for verification, including tni-miR- 



106a-3p, tni-miR-124b-3p, tni-miR-132-3p, tni-miR-142-3p, tni- 
miR-145b, tni-miR-17, tni-miR-223-5p, tni-miR-29d-3p, tni-miR- 
346, and tni-miR-378. The qPCR results showed the same up- 
and down-regulation patterns indicated by the sequencing data. 
Previous studies concerning the identified differentially expressed 
miRNAs showed that these miRNAs are involved in many 
immunological progresses in other species or cell lines. 

First, we focused on previous studies of miR-106a-3p, miR-17 
and miR-124b-3p in mammals. MiR-106a-3p, which belongs to 
the miR-17 family, is involved in many types of tumors [62,63,64]; 
additionally, miR-106a directly regulates Stat3 [65] and IL-10 
[66]. Stat3 is the signal transducer for IFN-y in some cell types 
[23], whereas IL-10 inhibits IFN-y production through inhibition 
of IL-12 [67,68]. MiR-17 is a member in the miR-17~92 cluster 
and acts as a potential oncogene [69] . In addition to its oncomiR 
function [70,71], miR-17 regulates macrophage differentiation 
[72] and macrophage inflammatory responses [73] and is also 
involved in Stat3 signaling [74]. Both miR-106a and miR-17 are 
regulated by c-Myc [75,76], which is a transcription factor that 
can be suppressed by IFN-y [77]. MiR-124b-3p belongs to the 
miR-124 family, and tni-miR-124b-3p is highly conserved with 
human miR-124 (100% identity to the mature sequence), which 
has been demonstrated to target Stat3 in different types of cancer 
[78,79,80]. These studies indicate that miR-106a-3p, miR-17, and 
miR-124b-3p may be involved in Stat3-dependent IFN-y signal- 
ing. MiR- 1 06a-3p may also be involved in feedback control of 
IFN-y expression by targeting IL-10, and miR-17 may also 
participate in the IFN-y-mediated activation of macrophages. 

We next focused on miRNAs relevant to IFN-y-induced 
cytokines. IL-12 is induced by IFN-y and then up-regxilates IFN- 
y through a feedback mechanism [81]. It has been reported that 
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Figure 6. Quantitative real-time PCR validation of selected miRNAs. (A): expression of miR-124b-3p, miR-142-3p, miR-145b, miR-17, miR-29d- 
3p, and miR-346 following IFNyl treatment, compared to the control group. (B): expression of miR-106a-3p, miR-124b-3p, miR-132-3p, miR-1 7, mlR- 
223-5p, and miR-378 following IFNy2 treatment, compared to the control group. The amount of each miRNA was normalized to that of U6 snRNA and 
is presented as the relative fold change (n = 3, mean ± SEM). Significant differences between the control and treatment groups are indicated (*P< 
0.05; **P<0.01; ♦**P<0.001). 
doi:10.1371/journal.pone.0096336.g006 



miR-1 32-3p is up-regulated after IL-12 treatment, and this up- 
regulation decreases Stat4 and IFN-y expression [82]. MiR-142- 
3p, a hematopoietic-specific miRNA [83], was down-regulated 
after rIFN-yl treatment in our study, and previous studies have 
shown that down-regulation of miR-1 42-3p can increase IL-6 
expression [84], resulting in CD4'*' T cell activation [85] and 
macrophage differentiation in tumors [86]. IL-6 is a proinflam- 
matory cytokine whose production can be activated by the LPS- 
induced macrophage response, and IFN-y is known to prime 
macrophages to enhance their response to LPS [87,88]. Addi- 
tionally, IFN-y primes helper T cell differentiation and macro- 
phage differentiation [89]. The suppression of miR-1 42-3p 
indicates a synergistic effect of miR-1 42-3p and IFN-y in the 
priming of macrophage responses. MiR-223-5p is the passenger 
strand of miR-223; although few studies about the passenger 
strand have been published, a decrease in miR-223 inhibits LPS- 
induced IFN-y production in splenic lymphocytes [90]. Addition- 
ally, down-regulation of miR-223 resulted in the activation of 
StatS and promoted the production of Toll-like-receptor (TLR)- 
triggered IL-6 and IL-lfi [91], indicating that miR-223-5p may 
function in the IFN-y priming of the LPS response and Stat3- 
dependent IFN-y signaling. Finally, although only a few studies of 



miR-346 have been undertaken, some reports have linked miR- 
346 to the immune system because of its ability to control the 
release of IL-18 and TNFot, and its targeting of TAP 1 [92,93,94]. 
These studies suggest that miR-346 may also be involved in the 
priming of the LPS response. 

Another down-regulated miRNA, tni-miR- 1 45b, is highly 
similar to human miR-1 45 (94.7% identical with one mismatch 
in the 3' end). Studies have shown that miR-1 45 expression is also 
decreased in various human tumors [95,96]. Additionally, in 
cancer cells, miR-1 45 is reported to target Statl [47], the main 
signal transducer of the IFN-y pathway, in addition to c-Myc 
[46,97]. Recent studies have also revealed that miR-1 45 regulates 
IFN-P expression by targeting the suppressor of cytokine signaling 
(SOCS)7 [98]. These data indicate that miR-145 may regulate tiie 
major IFN-y signaling pathway and may also participate in 
crosstalk among other interferons. Another miRNA, tni-miR-29d- 
3p, belongs to the miR-29 family and has 95.4% identity (one 
mismatch in the 3' end) with human miR-29b. The miR-29 family 
has been shown to direcdy target IFN-y [44,99], and IFN-y- 
induced Statl could control miR-29 expression through a 
feedback mechanism [45], suggesting that miR-29d-3p may be a 
crucial factor in the IFN-y-mediated immune responses that 
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Figure 7. KOG classes of predicted targets for Tetraodon nigroviridis miRNAs. 
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control the production of IFN-y itself. Previous studies have also 
shown that miR-378 can target the IFN-y receptor gene IFNGRl 
and suppress bovine luteal cell apoptosis [48]. These findings 
indicate that these differentially expressed miRNAs may play 
crucial roles in the IFN-y-mediated immune response, and the 
interactions between these miRNAs and IFN-y should be 
researched further. 

Additionally, seven out of the 10 vahdated miRNAs (miR-124b- 
3p, miR-132-3p, miR-142-3p, miR-145b, miR-17, miR-29d-3p, 
and miR-346) were difiFerentially expressed between the IFN-y 1 
and IFN-y2 treatment groups, suggesting functional segregation or 
synergy of the two fish IFN-y isoforms. Previous studies on the two 
fish IFN-y isoforms showed that they bind to different receptors 
and exhibit different expression patterns [36,40,41]. Comparative 
studies on the two IFN-y isoforms showed that they have 
complementary and also overlapping functions [39,40,41]. The 
miRNA differential expression data also suggest that the two 
isoforms have separate functions; for example, miR-124b-3p and 
miR-17 were down-regulated after rIFN-yl treatment but were 
significantly up-regulated after rIFN-y2 treatment. As mentioned 
above, miR-124b-3p and miR-17 may be involved in Stat3- 
dependent IFN-y signaling, which suggests that IFN-y 1 and IFN- 
y2 may have opposite effects on the Stat3 signaling pathway in T. 
nigroviridis. Previous studies in mammals have indicated that miR- 
1 06a and miR- 1 7 were up-regulated after 1 2 h of IFN-y treatment 
[100], which is in agreement with our findings for IFN-y2, 



suggesting that T. nigroviridis IFN-y2 may be closely related to 
mammalian IFN-y. Additionally, the rIFN-y 1 -specific down- 
re,gulation of miR- 145b indicates that IFN-y 1 may be involved 
in the Stat 1 -dependent signaling, and the rIFN-y 1 -specific up- 
regulation of miR-29d-3p sug,gests that IFN-y 1 may participate in 
the feedback control of IFN-y. Furthermore, the rIFN-y2-specific 
induction of miR- 1 32-3p indicates that IFN-y2 may be responsible 
for IL-12 production. The differential expression of miRNAs 
between the two treatment groups may provide important 
information regarding the functional segregation of the two fish 
type II IFN isoforms. 

The targets of all the T. nigroviridis miRNAs were predicted using 
Mireap software, and KOG analysis was performed. The KOG 
results showed that most of the target genes function in signal 
transduction mechanisms; translation, ribosomal structure, and 
biogenesis; and the cytoskeleton. Most of the target genes that 
function in signal transduction were also found in previous 
miRNA high-throughput sequencing studies [15,56]. To deter- 
mine the function of the differentially expressed miRNAs, their 
targets were predicted. Because the results of the target prediction 
were redundant, we chose only the top three targets with the 
highest scores for further analysis (Supplemental Table S7). These 
filtered targets were analyzed using GO enrichment and KEGG 
pathway analysis, and from the GO and KEGG analyses, we 
assembled a model of the miRNA-protein interaction network; 
however, the actual system may be much more complicated and 
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may require further study for elucidation. Some target genes were 
predicted to be involved in components of the IFN-y-mediated 
immune response, such as the Jak-Stat signahng pathway, the 
Toll-like receptor signaling pathway, the NOD-like receptor 
signaling pathway, the T cell receptor signaling pathway, antigen 
processing and presentation, and natural killer cell-mediated 
cytotoxicity. Targets involved in these pathways should be 
evaluated in future studies concerning IFN-y-induced miRNAs. 
However, computational methods for miRNA target prediction 
often have false positive rates [101,102]; thus, further analysis is 
needed to confirm the interaction between miRNAs and their 
targets. 

In conclusion, using high-throughput sequencing technology, 

we studied the T. nigroviridis miRNAs expressed in spleen cells. 
Many miRNAs were identified after rlFN-yl and rIFN-Y2 
stimulation of spleen c:ells, and dilferentially expressed miRNAs 
were identified and validated using qRT-PCR. The miRNA 
targets were predicted and analyzed using GO and KEGG 
analyses. This study provides the basis for the interaction between 
T. nigroviridis IFN-y and miRNAs, and these data provide clues for 
further studies on the mechanisms of IFN-y-mediated immune 
responses in T. nigroviridis. 
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